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Abstract: Fermi LAT has discovered two extended gamma-ray bubbles above and below the galactic plane. We propose 
that their origin is due to the energy release in the Galactic center (GC) as a result of quasi-periodic star accretion onto 
the central black hole. Shocks generated by these processes propagate into the Galactic halo and accelerate particles 
there. We show that electrons accelerated up to ~ 10 TeV may be responsible for the observed gamma-ray emission of 
the bubbles as a result of inverse Compton (IC) scattering on the relic photons. We also suggest that the Bubble could 
generate the flux of CR protons at energies > 10^^ eV because the shocks in the Bubble have much larger length scales 
and longer lifetimes in comparison with those in SNRs. This may explain the the CR spectrum above the knee. 
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1 Introduction 

Fermi bubbles are symmetric structures elongated above 
and below the Galactic plane for about 8 kpc 1 16 |. Their 
discovery is one of the most remarkable events in astro- 
physics. The origin of the bubble is still enigmatic and up 
to now a few models were presented in the literature. The 
team, which subtracted this structured gamma-ray emission 
from the total diffuse galactic emission presented different 
explanations of the phenomenon though they seem to trend 
towards the model of a single huge energy release in the 
Galactic center (GC) when about 10 million years ago a 
huge cloud of gas or a star cluster was captured by the cen- 
tral black hole that produced the Fermi bubbles seen today. 
Gamma-rays in their model are produced by the inverse 
Compton scattering of the electrons on relic photons. In 
IIT5I it was assumed that these electrons were accelerated 
in the bubble interior by a hypothetic MHD turbulence ex- 
cited behind the termination shock generated at the edges 
of the bubbles. Another interpretation was presented by 
0^ who proposed a relatively slow energy release 10'^^ 
erg s^^) due to supernova explosions as a source of proton 
production in the GC which emitted gamma-rays there. 
An alternative explanation based on the assumption that the 
energy for the Fermi bubbles was originated from star cap- 
ture events which occurred in the GC every 10"* — 10^ years 
was suggested by [9 |. About one thousands of such events 
form giant shocks propagating through the central part of 



the Galactic halo and thus produce accelerated particle re- 
sponsible for the bubble emission. 

Processes of particle acceleration by the bubble shocks in 
terms of sizes of the envelope, maximum energy of accel- 
erated particles, etc. may differ significantly from those 
obtained for SNs that may lead to the maximum energy 
of accelerated particles much larger than can be reached 
in SNRs. In this respect, we assume that acceleration of 
protons in Fermi bubbles may contribute to the total flux 
of the Galactic cosmic rays (CR) above the "knee" break 
(> IQi^eV). 

2 Structure of Shocks in the Fermi Bubble 

As it was assumed in fO) the central massive black hole 
captures a star every tq ~ 10'* — 10''' years and as a re- 
sult releases energy about £q ^ 10^^ erg in the form of 
subrelativistic particles which heat the central < 100 pc in 
the GC. This heating produces a shock propagating into the 
surrounding medium [18 |. For an exponential atmosphere 
with the plasma density profile 



p{z) = po exp 
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an analytical solution for permanent energy injection by 
star accretion was obtained by lfT3l . Here z is the coor- 
dinate perpendicular to the Galactic plane. The radius of 
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Figure 1; The bubble multi-shock structure, 
the bubble as a function of the height z and the time t is 
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y is a current volume enveloped by the shock 

a{t) 

V{t) = 27r / r^{z,t)dz, 
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a is the position of the shock top 
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W — f o/tq is the average luminosity of the central source, 
7g is the polytropic coefficient, and a and A are numbers 

m. 

Depending on capture parameters one can imagine the bub- 
ble interior as a volume fiUed with many shocks of different 
ages (see Fig[l])- 

3 Particle Acceleration by the Bubble 
Shocks 

Correct analysis of shock acceleration in the bubble re- 
quires sophisticated calculations of each stage of this pro- 
cess which we hope to perform later. Now we present 
simple estimates of characteristics of accelerated spectra. 
We analyze spectra of protons and electrons separately be- 
cause they are formed by completely different processes 
and as we later conclude their spectral characteristics dif- 
fer strongly from each other. We present also the spectrum 
of gamma-ray emission produced by the inverse Compton 
(IC) scattering of the accelerated electrons in the bubble. 



3.1 Proton Spectrum 

Below we analyse the spectrum of protons accelerated in 
the bubble and discuss whether the bubble contribution 
to the total flux of CRs in the Galaxy may explain the 
knee steepening. We remind that the generally accepted 
point of view is that the flux of relatively low energy CRs 
(< lO^^eV) is generated by SNRs which ejects a power- 
law spectrum E^'^ into the interstellar medium. This spec- 
trum is steepened by propagation (escape) processes in 
the Galaxy in accordance with the spectrum observed near 
Earth (for details see |6|). However these sources can 
hardly produce CRs with energies more than lO^^eV (see 
imiSj). Just at this energy a steepening (the knee) in the CR 
spectrum is observed. We assume that the bubble may gen- 
erate the flux of CRs at energies above lO^^eV. The origin 
of CRs with energies above lO^^eV is the goal of this sub- 
section. It is natural to assume that the bubble shocks may 
produce this flux considering thek large scales and long 
lifetimes. 

In the framework of our model one can imagine the bubble 
as at cylinder of the radius ~ 3 — 6 kpc filled with hundreds 
of shocks propagating in series one after another (see Fig. 
[U. The frequency of shock injection (star capture) is poorly 
known parameter, especially for the conditions of the GC. 
From the analysis of shock hydrodynamics we estimated 
the age of outer shock by 10^ yr though this value depends 
strongly on parameters of the self-similar solution which is 
an idealization of conditions in the Galaxy. The frequency 
of star capture by a central black hole is estimated from 
theoretical treatments v ^ 10^* — 10~^yr~^ |i2|. Thus, a 
multi-shock structure can be produced by the capture pro- 
cesses with an average distance L between separate shocks 
given by 

L = TcapU = iO{v X 3 X I0^yr){u/\{fcm/s)pc. (6) 

On the other hand there is another spatial scale which 
characterizes processes of particle acceleration by a single 
shock which is Id ^ D/u where D is the spatial diffu- 
sion coefficient near a shock and u is the shock velocity. 
Depending on the relation between L and Id there may 
be different regimes of acceleration in the bubble. This 
problem of particle acceleration in supersonic turbulence 
or multi-shock structure was analyzed in fSl. The accelera- 
tion regime is characterized by a dimensionless parameter 
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If "0 << 1, the diffusion length scale of a single shock 
Id exceeds the distance between shocks L. Therefore, 
particle are accelerated by interactions with many shocks 
that gives the acceleration similar to the classical stochastic 
Fermi acceleration. The critical energy Ei for this regime 
of acceleration is estimated from the condition ij} ^ \ or 
I (El) ^ L. In the Bohm limit with the diffusion coeffi- 
cient D ^ ctl (E) 1 3 where (E) = E/eB is the particle 



32nd International Cosmic Ray Conference, Beijing 201 1 



Larmor, the value of Ei is 



El 



eBLu 



10 



15 



lO^cm/s 



eV. 
(8) 

which is about the position of the knee break. 

The kinetic equation for the case tp << 1{E >> Ei) only 
contains the diffusion in the spatial and momentum spaces 
(Fermi stochastic acceleration) and has the form 
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(9) 



where p and z are the cylindrical spatial coordinates, p is 
the particle momentum. D{p,p) is the spatial diffusion co- 
efficient and k{p, p) is the momentum diffusion coefficient. 



(10) 



The term Q{p,z,p) describes CR injection by supernova 
remnants in the disk with energies E < 10^^ eV and with 
observed power-law spectrum. 

The momentum dependence of / can be presented by a 
power-law function, f{p) oc p^'^, where 7 should be de- 
termined from Eq.(|9]i: 
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where H is the height of the Galactic halo and D = 
Da{E/Eo)^ is the average diffusion coefficient in the 
Galaxy. For ~ 10 kpc, Dqc^Sx lO^^ cm^s"^ ^ ~ 0.3, 
~ 3 GeV (n] and w 10^ cm/s, L = 30 pc we have 
7 ~ 5, i.e., the power-law spectrum above 10^^ eV is in 
the form N{E) oc E~^ as necessary for the knee spec- 
trum. The proton spectrum derived from Eq. (|9|l is shown 
in Fig. |2] 

3.2 Electron Spectrum and Spatial Distribution 

Unlike protons relativistic electrons lose their energy effec- 
tively by synchrotron and inverse Compton radiation. The 
rate of energy loss of electrons is 



dE 
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The maximum energy of electrons E^ax can be estimated 
if the spatial diffusion coefficient at the shock Dgh is 
known. From Eq. (fT2l i we have 
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For the Bohm diffusion coefficient we can estimate the up- 
per limit for E'^-^^^ which is 
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Figure 2: The Bubble contribution to the flux of CRs. Ti- 
bet, KASCADE, HiRes and Auger data are summarized in 
[14|. ATIC data are from [IJ, Proton data are from fT2l . 
RUNJOB data are from iH. 



that gives for the velocity u = 10* cm s~^, the magnetic 
field strength H = 10"*^ G and Wph = 0.25 eV cm"^ the 
maximum energy of accelerated electrons about EBmax ~ 
5 X lO^^eV, i.e., for electrons E < Esmax « Ei. It 
follows from this inequality that electrons are accelerated 
in the regime ^ >> 1 (single shock acceleration). Kinetic 
equations for this case were derived in |[8l. 

In the case of electrons there is one more spatial parameter 
essential for acceleration by shocks, namely the electron 
mean free path A. Therefore we have two more dimension- 
less parameters which describe electron acceleration in the 
case of supersonic turbulence, Ish/^ and Id/^- Nearby 
shocks electrons propagate by diffusion, and A there is a 
function of energy 



XniE) 




(15) 



From Eq. ( fT3] l one can see that for E < Ef^^^.^ we have 
A > Id, i.e., acceleration by shocks is effective. 

Far away from shocks electrons propagate by convection 
and for the velocity u the mean free path is 



XviE) 



I3E 



(16) 



For relatively low energy electrons X{E) > Igh, and a 
regime of multi-shock acceleration for electrons is realized 
in this case where interactions with many shocks change 
slowly the shock to E~^ (for details see lH)) which pro- 
duces a hard E~^ spectrum . The break position E* be- 
tween the single shock and multi-shock spectra follows 
from the equality \{E*) = Igh 



E* 



pish 



(17) 



that gives, e.g., E* - 2.8 x 10" eV for Uh = 10^^ cm. 
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If we assume that the bubble gamma-rays are produced by 
IC scattering of electrons on the relic photons, then it fol- 
lows from our model that: 1 . a drop of the bubble gamma- 
ray flux ai < 1 GeV as observed by IJ6i is due to a 
flattening of electron spectrum at _E < E*; and 2. a drop 
of this spectrum in the range E^ > 100 GeV is due to a 
cut-off in the electron spectrum at E ^ ^max- The ex- 
pected spectrum of gamma-ray emission from the Bubble 
due to IC scattering of the electrons is shown in Fig. |3] 



10- 




i 1-4- 




10 



10° 



10' 



10'" 
E(eV) 



10' 



10" 



Figure 3: The spectrum of gamma-ray emission from 
Fermi bubble in case of multi-shock acceleration. Data 
points are taken from 1 16 1. 
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Figure 4: Spatial distribution of the gamma-ray emission. 
Data are from fT6l . Top: in case of single shock. Dotted 
line correspond to D=10^^ cm^/s, solid to D=10'^'^ cm^/s, 
dashed to D=10'^^ cm^/s and dash-dotted D=10''^ cm^/s, 
Bottom: in case of several shocks distributed in accordance 
with 

In Fig. |4] we show an expected spatial distributions of 
gamma-ray emission from the bubble for the single shock 
and multiple shocks cases. For the multiple shocks it was 
assumed that they distributed in accordance to the solution 
(|2|l and it was also assumed that the amount of accelerated 
electrons is proportional to the density of shocks. From this 
figure one can see that the single shock model (without the 
second order Fermi acceleration in the bubble interior) is 



unable to reproduce the data. However, for the parameters 
of star capture model these data are nicely described. 

4 Conclusion 

We have shown that series of shocks produced by a se- 
quential stellar captures by the central black hole can fur- 
ther re-accelerate the protons emitted by SNRs up to en- 
ergies above lO^^eV. The predicted CR spectrum con- 
tributed by the Bubble may be i?^" where i' ^ 3 for 
10^^ eW < E < 10^^ eV that explains the knee CR spec- 
trum. 

The regime of electron acceleration in the bubble is quite 
different from that of protons. It is a combination of single 
and multishock accelerations. In this case we have a cut-off 
of the electron spectrum at > 3 10^'^ eV and flattening 
of the spectrum at E < 100 GeV that explains nicely the 
bubble gamma-ray spectrum and the sharp edge spatial dis- 
tribution observed by 1,16,1 if this emission is due to IC on 
the relic photons 
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